Spatial persistence and survival probabilities for fluctuating interfaces 
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We report the results of numerical investigations of the steady-state (SS) and finite-initial- 
conditions (FIC) spatial persistence and survival probabilities for (H-l)-dimensional interfaces with 
dynamics governed by the nonlinear Kardar-Parisi-Zhang (KPZ) equation and the linear Edwards- 
Wilkinson (EW) equation with both white (uncorrelated) and colored (spatially correlated) noise. 
fH ^ We study the effects of a finite sampling distance on the measured spatial persistence probability and 

^ . show that both SS and FIC persistence probabilities exhibit simple scaling behavior as a function 

O ' of the system size and the sampling distance. Analytical expressions for the exponents associated 

with the power-law decay of SS and FIC spatial persistence probabilities of the EW equation with 
power-law correlated noise are established and numerically verified. 
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I. INTRODUCTION 
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Q ■ The concept of temporal persistence which is closely related to first-passage statistics, has been used recently 
' to study various non-Markovian stochastic processes both theoretically 2, 3] and experimentally Another 
quantity of interest in the study of the statistics of spatially extended systems is its natural analog, the spatial 
' persistence probability. This idea has been investigated theoretically in the context of {d + 1) -dimensional Gaussian 
^ interfaces with dynamics described by linear Langevin equations, where the variable undergoing stochastic evolution 
T— I is the height h{x, t) of the interfacial sites {x is the lateral position along the interface and t is the time). The spatial 
\^ [ persistence probability of fluctuating interfaces, denoted by P{xq,xo -\- x), is simply the probability that the height 
\l • of a steady-state interface configuration, measured at a fixed time Iq, does not return to its "original" value /i(xo,io) 
I at the initial point xq within a distance x measured from xq along the interface. In the long-time, steady-state limit, 
. the spatial persistence probability P{xq,xo + x), which depends only on x for a translationally invariant interface, 
' has been shown Q to exhibit a power-law decay, P(a;o,a;o -I- a;) ~ x~^ . One of the interesting results reported in 
, Ref. is that the spatial persistence exponent 9 can take two values determined by the initial conditions or selection 
rules imposed on the starting point Xq: 1) 9 = 9ss, the "steady state" (SS) persistence exponent if xq is sampled 
^ , uniformly from all the sites of a steady-state configuration; and 2) 9 = 9pic^ the so-called finite-initial- conditions 
I ■ (FIC) persistence exponent if the sampling of xo is performed from a subset of steady-state sites where the height 
' O ' variable and its spatial derivatives are finite. The spatial persistence probabilities obtained for these two different 
^ . ways of sampling the initial point are denoted by Pss{xo, xq -\- x) and Pfic{xq,xq -\- x), respectively. 

The values of the exponents 9ss and 9fic for interfaces with dynamics described by a class of linear Langevin 
equations have been determined in Ref. P| using a mapping between the spatial statistical properties of the interface 
in the steady state and the temporal properties of stochastic processes described by a generalized random-walk 
equation. It turns out that for these systems, 9ss is equal to either 3/2 — n for 1/2 < n < 3/2 or for n > 3/2, 
5^ • where n = {z — d -\- l)/2, d is the spatial dimension, and z is the standard dynamical exponent of the underlying 
' Langevin equation. The FIC spatial persistence exponent is found to have the value 9fic — where 9{n) is a 

temporal persistence exponent for the generalized random walk problem to which the spatial statistics of the interface 
is mapped. Two exact results for 9{n) are available in the literature: 9{n = 1) = 1/2, corresponding to the classical 
Brownian motion 9] and 9{n = 2) = 1/4, corresponding to the random acceleration problem ,10|. 

Very recently, experimental measurements of the spatial persistence probability have been performed for a 
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system (combustion fronts in paper) that is believed to belong to the Kardar-Parisi-Zhang (KPZ) [Tl| universality 
class. However, the FIC spatial persistence probability is not investigated at all in this work. Instead, the authors 
analyze a "transient" spatial persistence (i.e., the probability is measured by sampling over all the sites of a transient 
interfacial profile obtained before the steady state is reached). This transient spatial persistence is completely different 
from the FIC spatial persistence which is measured in the steady-state regime by sampling a special class of initial sites. 
As a consequence, additional study is required in order to understand the experimental and numerical possibilities 
for measuring Ppic and its associated nontrivial exponent 9fic- 

In this paper, we present the results of a detailed numerical study of spatial persistence in a class of one-dimensional 
models of fluctuating interfaces. Our interest in analyzing the spatial persistence of fluctuating interfaces is motivated 
to a large extent by their important (and far from completely understood) role in the rapidly developing field of 
nanotechnology where the desired stability of nanodevices requires understanding and controlling thermal interfacial 
fluctuations. In this context, the study of first-passage statistics in general, or of the persistence probability (both 
spatial and temporal) |S,|3J in particular, turns out to be a very useful approach. To address this problem we consider 
stochastic interfaces with dynamics governed by the Edwards-Wilkinson (EW) and KPZ equations. For the EW 
equation, we consider both white noise (uncorrelated in both space and time) and "colored" noise that is correlated in 
space but uncorrelated in time. The effect of noise in spatially distributed systems is an interesting problem by itself 
and has been widely studied In this paper, we investigate the effects of noise statistics on the spatial structure 
of fluctuating interfaces using the conceptual tool of spatial persistence probability. Using the isomorphic mapping 
procedure of Ref. we derive exact analytical results for the spatial persistence exponents of {d + 1) -dimensional 
EW interfaces driven by power-law correlated noise. We then compare our analytical results with those obtained from 
numerical integrations of the corresponding stochastic equations. The use of power-law correlated noise in the EW 
equation allows us to explore the situation where the two spatial persistence exponents 9s s and Opic are different. 

Our numerical study also provide a characterization of the scaling behavior of spatial persistence probabilities as 
functions of the system size. Information about the system-size dependence of persistence probabilities is necessary 
for extracting the persistence exponents from experimental and numerical data. In studies of the scaling behavior of 
spatial persistence probabilities, one has to consider another important length scale that always appears in practical 
measurements: this is the sampling distance Sx which represents the "nearest-neighbor" spacing of the uniform grid 
of spatial points where the height variable /i(a;,io) is measured at a fixed time Iq. The sampling distance Sx is the 
spatial analog of the "sampling time" [l4l ITsI l that represents the time-interval between two successive measurements 
of the height at a fixed position in experimental and computational studies of temporal persistence. Once the effect of 
a finite 6x on the measured spatial persistence is understood, one can relate correctly the experimental and numerical 
results to the theoretical predictions. Our study shows that the spatial persistence probabilities (both SS and FIC) 
exhibit simple scaling behavior as functions of the system size and the sampling distance. 

In addition to the temporal persistence probability, the temporal survival probability [sl llSj has been shown recently 
to represent an alternative valuable statistical tool for investigations of first-passage properties of spatially extended 
systems with stochastic evolution. In the context of interface dynamics, the temporal survival probability is defined as 
the probability that the height of the interface at a fixed position does not cross its time-averaged value over time t. In 
contrast to the power-law behavior of the temporal persistence probability (which, we recall, measures the probability 
of not returning to the initial position), the temporal survival probability exhibits an exponential decay at long times, 
providing information about the underlying physical mechanisms and their associated time scales I5j. In this study, 
we make the first attempt to analyze the behavior of the spatial survival probability, S{xo,xo + x), defined as the 
probability of the interface height between points xq (which is an arbitrarily chosen initial position) and xq + x not 
reaching the average level (h) (rather than the original value h{xo, to)). We present numerical results for S(xo, xq -\- x) 
that show that its spatial behavior in the SS regime is neither power-law, nor exponential, while in the FIC regime, 
it becomes very similar to the spatial persistence probability, Pfic{xo, Xq + x). 

The paper is organized as follows. In Sec. m we define the models studied in this paper, review existing analytical 
results about their spatial persistence properties, and present new analytical expressions for the spatial persistence 
exponents for EW interfaces with colored noise in arbitrary spatial dimension. In Sec. IIIII we describe the numerical 
methods used in our study and discuss how the spatial persistence and survival probabilities are measured in our 
numerical simulations. The results obtained in our (l-l-I)-dimensional numerical investigations are described in detail 
and discussed in Sec. IIVI for both discrete stochastic solid-on-solid models (Sec. llVt V) and the spatially discretized 
EW equation with colored noise (Sec. IIVB ). Sec. IVl contains a summary of the main results and a few concluding 
remarks. 
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II. STOCHASTIC EQUATIONS FOR FLUCTUATING INTERFACES 



We have performed a detailed numerical study of the spatial persistence of (l+l)-dimensional fluctuating interfaces 
where the dynamics is described by the well known EW equation 

^!^ = V^hix,t)+v{x,t), (1) 

or alternatively by the KPZ equation 

dh{x,t) 



dt 



Vh{x,t) + {\7h{x,t)y +Tj{x,t), (2) 



where V and refer to spatial derivatives with respect to x, and r]{x,t) with {ri{x,t)ri(x ,< )) cx 5{x — x )6{t — t ) 
is the usual uncorrelated random Gaussian noise. The dynamical exponent for Eq.(^ is z = 2, and since d = 1 in 
our study, the variable n defined in Sec. QJis equal to 1. So, we expect both 9ss and 9fic for this system to be 
equal to 1/2 0. Although the KPZ equation is nonlinear, characterized by z = 3/2, it is well-known that in the long 
time limit, the probability distribution of the stochastic height variable h{x, t) in this equation is the same as that in 
the EW equation (i.e. P{h) ^ exp[ — J dxiyhY ] ) in (1+1) dimensions. The static roughness exponent, a, is the 
same (a — 1/2) for both cases. The 1+1-dimensional KPZ model differs from the EW model in the transient scaling 
regime where the interfacial roughness grows as a power-law in time, but this temporal regime is not involved in the 
calculation of the spatial persistence probabilities, as explained in Sec.^ As a consequence, the steady-state spatial 
properties of (l-fl) -dimensional interfaces governed by Eq. ((SJ can be mapped, as for Eq. into a stochastic process 
with n=\. So, the expected values oiOss and 9fic for the (l-l-l)-dimensional KPZ universality class are also equal 
to 1/2. Thus, studies of (1-1-1) -dimensional KPZ and EW interfaces do not bring out the interesting possibility of 
different values for the spatial persistence exponents Oss and Opic- 

To examine the theoretical prediction Q of a possible difference between the values of Oss and Opic, we consider 
the case when the interface dynamics is governed by a EW-type equation with long-range spatial correlations in the 
noise. Specifically, we consider Eq. ^ with Gaussian colored noise [T^ with variance given by 

{'nc{x,t) Ti,{x',t))^gp{x-x')5{t~t'), (3) 
where < P < 1/2 is a parameter that characterizes the spatial correlation of the noise, and 

, „ r \x-x'\^P-^ if \x-x'\^0 

We have chosen gp{0) as in Ref. |0| (i.e. gp{0) — l//9(l/2)^^). As discussed below, the SS and FIG spatial persistence 
exponents for (l-l-l)-dimensional interfaces described by the EW equation with this kind of colored noise are expected 
to be different from one another. This system, thus, provides an opportunity to examine in detail the role of the 
choice of the initial points in determining the form of the decay of the spatial persistence probability. 

By applying the isomorphic mapping recipe of Ref. to the {d + l)-dimensional version of Eq. with colored 
noise r]c whose statistics is defined by Eqs. (j3Jl and Q), we obtain the result n = (z — d+l)/2 + p with z = 2, implying 
the following analytical expressions for the spatial persistence exponents: 

Oss^^-P (5) 



and 



Ofic = 0['^+p]- (6) 



Thus, the value of 9ss is completely determined by the noise correlation parameter p. However, based on the range 
of values for p, we can only infer that Opjc varies (presumably in a continuous manner) between 9{^^) and ^(^-^) 
as the parameter p is increased from to 1/2. For d — 1, this implies a change from the value 9{1) — 1/2 to 0(3/2), 
expected to lie between 1/2 and 9{2) = 1/4, as p changes from to 1/2. Since the value of 9ss for d = 1 goes to as 
p approaches the value 1/2, it is clear that the values of the two spatial persistence exponents must be different for a 
general value of p in the range [0,1/2). This difference would be small for p near zero (the two persistence exponents 
have the same value for p = 0), and maximum for p near 1/2. Therefore, the model with p substantially different 
from zero provides a numerically tractable situation where the interesting theoretical prediction of the existence of 
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two different nontrivial spatial persistence exponents can be tested. We also mention that the usual dynamical scaling 
exponents take the following p-dependent values in the model with colored noise: a — (2 — c?+2p)/2, /3 = (2 — d+2p)/4. 
Thus, the general result 0, 9ss = 1 — a, is satisfied for all d and p. 

We have investigated these aspects in a detailed numerical study of models that belong in the universality classes 
of the Langevin equations of Eqs. and ((SJ- For Eq. (QJ with uncorrelated white noise, we have used a discrete 
stochastic solid-on-solid model (the Family model [l3|) which is rigorously known to belong to the same dynamical 
universality class. For Eq. |(2Jl with uncorrelated white noise, we have also used a discrete solid-on-solid model (the 
Kim-Kosterlitz model 18]). Finally, for the EW equation with colored noise, the numerical results were obtained 
from a direct numerical integration of the spatially discretized stochastic differential equation. 



III. NUMERICAL METHODS 



Simulations of the atomistic Family and Kim-Kostcrlitz models are carried out using the standard Monte Carlo 
method for implementing the stochastic deposition rules of each model. Numerical integration of the EW equation 
with colored noise is performed using the simple Euler method We solve the (l-)-l)-dimensional Eq. (1} with 

spatially long-range correlated noise for the real variable h{xj,tn), where t„ — nAt {n = 0, 1, . . .) and Xj — jAx 
(j = 0, 1, . . . , L — 1), with periodic boundary conditions. Here, At and Ax are the spatial and temporal grid spacings, 
respectively. Using the forward-time centered-space representation 19], Eq. becomes: 



h{Xj,tn+l) - h{Xj,ta) = At 



h{Xj+i,tn) - 2h{Xj,tn) + h{Xj-i,t„) 

{Axy 



+ VAtric{xj,tn). (7) 



We have chosen Aa; — 1 and At small enough (i.e. At — 0.01) in order to satisfy the stability criterion 2At/ (Ax)^ < 1. 
The spatial correlation of the noise is given by 

{T]c{Xj,tn) r]cixk,tm)) =gp{Xj - Xk)Sn,Tn (8) 

with 

( \xj-Xk\'^P-^ if l<\xj-Xk\<^ 
9p{xj - Xk) = < (i - \xj - Xk\fP-^ if \xj - Xfc] > f (9) 
[ 5p(0) if Xj - Xfe = 

where gp{Q) = l/p{l/2)'^P. The colored noise is generated using the recipe from Ref. The Fast Fourier Transform 
operation that is used in the noise- generation procedure constrains the system size to be an integral power of 2. Due 
to the use of periodic boundary conditions (which are also imposed on the noise correlation function, see Eq. Q), 
the range of x over which spatial correlations and persistence properties are meaningfully measured is of the order of 
L/2. 

The SS spatial persistence probability Pssixoi xq + x) is measured at a fixed time to (which is much larger than the 
time tgat ~ required for the interface roughness to saturate) as the probability that the interface height variable 
does not cross its value, h{xo,to), at the initial point Xq as one moves along the interface from the point xq to the point 
Xq + X. This probability is averaged over all the sites in a steady-state configuration and also over many independent 
realizations of the stochastic evolution. Thus, 

Pss{xoiXq + x) = Prob { sign [h{xo + x') — h{xo)] = constant, VO<a;'<a;, VxqG Sss Ii (10) 

where sign [y] represents the sign of the fluctuating quantity y, and Sss is the ensemble containing all the lattice 
sites in a steady-state configuration. The FIC spatial persistence probability Ppic {xo ,xo + x) is obtained in a similar 
manner, except that the average is performed over a particular subensemble of the steady-state configuration sites, 
Sfic C SsSi characterized by finite values of the height variable and its spatial derivatives: 

PpicixQ, Xq + x) = Prob { sign [h{xQ + x') — h{xo)] = constant, \f < x' < x , VxqG Sfic }■ (H) 

Since the persistence probabilities are averaged over the choice of the initial point xq , we omit writing xq explicitly 
in the arguments of Pss and Ppic from now on, while stressing the important fact that the ensemble of initial sites 
used in the averaging process determines which one of the two persistence probabilities is obtained. We consider 
two different methods for measuring Ppic (x) , depending on the type of the model (atomistic solid-on-solid model or 
spatially discretized Langevin equation) being studied. In the former case where the height variables are integers, the 
FIC spatial persistence probability measurement involves a sampling procedure from the subset of sites characterized 
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by a fixed integer value of the height (measured from the average, (h), of the heights of all the sites at time to) which 
is substantially smaller than the typical value of the height fluctuations measured by the saturation width of the 
interface profile. In calculations using the direct numerical integration technique, the height variable can take any 
real value. So, the probability of finding a fixed value of the stochastic height variable is infinitesimally small. For 
this reason, fixing a reference level H and sampling over the sites with h{xo, tg) — (h) + H is useless. We, therefore, 
consider in this case a continuous interval of height values (symmetric with respect to the average height (h)) with 
width w which is considerably smaller than the amplitude of the height fluctuations. The positions characterized by a 
height variable within this interval represent the subensemble of lattice positions involved in the sampling procedure 
necessary for measuring Ppic{x). 

The spatial survival probabilities corresponding to the SS and FIC conditions are calculated similarly to the corre- 
sponding persistence probabilities, except that the stochastic variable under consideration becomes h[xQ + x') — (h). 
Thus, 

Sssi^o^ xq + x) = Prob { sign [h{xQ + x') — {h)] = constant , VO<x'<a;, M xq & Sss }i (12) 

and 

SFicixQ^xa + x) = Prob { sign [h(xQ + x') — {h)] — constant VO<a:'<a;, VxoG Spic }■ (13) 

IV. RESULTS AND DISCUSSIONS 
A. Solid-on-solid models 

In the solid-on-solid Family and Kim-Kosterlitz models, the interface configuration is characterized by a set of 
integer height variables {/ii}i=i,L corresponding to the lattice sites i = 1, . . . , L, with periodic boundary conditions. 
Since all the measurements of the spatial persistence and survival probabilities are done in the steady-state regime (i.e. 
in the regime where the interfacial roughness has reached a time-independent saturation value), we used relatively 
small systems with L ~ 200 — 3000 in order to be able to achieve the the steady state within reasonable simulation 
times. The resulting steady-state interfacial profile, corresponding to a final time to >> L^, is used to compute the 
spatial persistence and survival probabilities. The calculation of Pss{x) is relatively simple: it involves measuring the 
fraction of initial lattice positions (all possible choices of the initial point are allowed) for which the interface height 
has not returned to the height of the initial point (for persistence probability) or to the average height level (h) (for 
survival probability) over a distance x, averaged over many independent realizations 10^ — 10**) of the steady state 
configuration. Measurements of Ppici^) oi' Spjc{x) involve, in addition to these steps, a preliminary selection of a 
subensemble of lattice sites which are characterized by a fixed and small value H of the height measured relative to 
the spatial average. Only the sites that belong to this subensemble (i.e. only the sites with hi = H + (h)) are used 
as initial points in the FIC measurements. 

Two distinct length scales have to be taken into consideration in the interpretation of the numerical results for the 
spatial persistence probability: the size L of the sample used in the simulation, and the sampling distance Sx which 
denotes the spacing between two successive points where the height variables are measured in the calculation of the 
persistence probability. The minimum value of Sx is obviously one lattice spacing, but one can use a larger integral 
value of Sx in the calculation of persistence and survival probabilities. For example, a calculation of the persistence 
probability with Sx — m would correspond to checking the heights of only the sites with index io + jm, where io 
is the index of the initial site and j = 1,2, . . .. While the importance of L in the measurement of P{x) is obvious 
(it sets the maximum distance for which P{x) can be meaningfully measured), the effect of Sx is rather intricate 
and has to be carefully investigated. In Fig.^ we start to analyze these effects by looking at Pss{x) for EW-type 
interfaces. We note that when Pssi^) is measured in systems with different sizes, using the smallest possible value 
for dx (i.e. 6x = 1), the exponent associated with the power-law decay of the persistence probability does not change, 
but there is an abrupt downward departure from a power-law behavior near x = L/2. It is not difficult to understand 
this behavior qualitatively: as discussed earlier, measurements of spatial correlations and persistence probabilities 
in a finite system of size L with periodic boundary conditions are meaningful only for distances smaller than L/2. 
In Fig. we have shown the results for Pssi^) when L remains fixed and 6x is varied. Since the the persistence 
probability is, by definition, equal to unity for x = Sx (see Eq. (|10() '). we have plotted Pss as a function of x/Sx in 
this figure to ensure that the plots for different values of 6x coincide for small values of the x-coordinate. The plots 
for different dx are found to splay away from each other at large values of x/Sx, with the plots for larger Sx exhibiting 
more pronounced downward bending. Again, the reason for this behavior is qualitatively clear: since a double-log 
plot of Pss{x) vs X begins to deviate substantially from linearity as x approaches L/2 (see Fig. ^), the downward 
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FIG. 1: The steady state spatial persistence probability, Pssix), for (l+l)-dimensional EW interfaces with white noise, 
obtained using the discrete Family model. Panel (a): Double-log plots of Pssix) vs x for a fixed sampling distance Sx — 1, 
using three different values of L, as indicated in the legend. Panel (b): Double-log plots of Pss{x) vs x/Sx for a fixed system 
size, L = 1000, and three different values of Sx, as indicated in the legend. 

bending of the plots in Fig. ^ (which are all for a fixed value of L) occurs at a smaller value of x/Sx for larger Sx. 
A more detailed scaling analysis of the dependence of the persistence probabilities on x and Sx is described below. 

In Fig. 121 we show the results for spatial persistence and survival probabilities for the discrete Family model. It is 
obvious from the plots that the spatial persistence probabilities Pssix) (panel (a)) and Pfic{x) (panel (c)) exhibit 
power-law decays over an extended range of x values. The abrupt decay to zero near x — L/2 is due, as discussed 
above, to finite size effects. The spatial persistence exponents are extracted from the power-law fits shown in the 
log-log plots as dashed straight lines. We find that 6gg ~ 0.51, in good agreement with the expected value 1/2. 
However, it is clear that the steady state survival probability Sss{x), shown in Fig. does not exhibit a power-law 
behavior. This is similar to the qualitative behavior of the temporal survival probability in the steady state of the 
Family model p5l |. 

We now return to the dependence of the persistence probabilities on the sample size L and the sampling distance Sx. 
Since L and Sx are the only two length scales in the problem (the lattice parameter serves as the unit of length) , it is 
reasonable to expect that the persistence probabilities would be functions of the (dimensionless) scaling variables 
x/L and Sx/L. If this is true, then plots of P vs. x/L for different sample sizes should show a scaling collapse if the 
ratio Sx/L is kept constant. A similar scaling behavior of the temporal survival probability as functions of L and 
the sampling time St (in that case, the scaling variables are t/L^ and St/L^) was found in Ref. indicated in 

panels(b-d) of Fig. |21 we have used various values for the sampling distance Sx in the measurement of Pssi^) and 
Ppici^)- We observe that when the sampling distance is increased in proportion to the system size (so that Sx/L is 
held fixed), all the Pssi^) curves collapse when plotted vs. x/L (see panel (b)). This confirms that the scaling form 
of the steady state persistence probability is: 

Pssix,L,Sx)^ h{x/L, Sx/L), (14) 

where the function fi(xi,X2) shows a power-law decay with exponent 6ss as a function of xi for small values of xi 
and a;2 <C 1. 

Let us turn our attention to Ppic (x) . In the data shown in panel (c) of Fig. |21 we have chosen the subensemble 
Spic of sampling positions to contain only the lattice sites whose height hi is equal to the average value (h) (i.e. 
H = 0). Obviously, in this case the definitions for persistence and survival probabilities become identical, since the 
probability that the height variable does not return to the original value (i.e. hi = (h)) is precisely the probability 
that the height variable does not reach the average level (h). We find that 9fic — 0.48 using a system with L = 1000 
and Sx — 1 and considering the subensemble of sites with H — 0. We note that a remarkable collapse of Pfic{x) vs. 
x/L curves for different values of L is again obtained when Sx is adjusted to be proportional to the system size L, as 
shown in panel (c). More interestingly, we observe that fixing the level H to a nonzero value introduces a "height" 
scale in the problem that is related to the steady-state value of the interface width. Since this width is proportional 
to L", where a is the roughness exponent, we expect the dependence of Ppic on H for nonzero values of H to be 
described by the scaling variable H/L"-. We observe that if the level H is chosen to be proportional to L", then the 
calculated values of Ppic for different sample sizes, obtained using values of Sx such that the ratio Sx/L is also held 
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FIG. 2: The spatial persistence probabilities, Pss(x) and Pfic{x), and the spatial survival probability, Sssix), obtained from 
simulations of the Family model in (1+1) dimensions. In panels (a) and (b) we show the data for Pssix) and Sssix), while 
in panels (c) and (d) we display the data for Pfic{x). Panel (a): Pss{x) and Sssix) for L = 1000, 5x — 1. The dashed 
line represents the best fit of the Pssix) data to a power-law form. Panel (b): Finite-size scaling of Pss{x, L, 5x). Three 
probability curves are obtained for three different sample sizes with the same value for the ratio Sx/L = 1/200. Panel (c): 
Scaling of Pfic{x, L, Sx, H) for the same values of L and Sx as in panel (b). Pfic is calculated by sampling over lattice sites 
with H = 0. Panel (d): Scaling of Pfic{x, L, 5x, H) for three different sample sizes with the same value for the ratio Sx/L, 
sampling over two subsets of lattice sites with the same value of H/L°' {a = 0.5): l/-\/200 (upper plot) and 4/\/200 (lower 
plot). 



constant, exhibit a perfect scaling collapse, as shown in panel (d) of Fig.|21 This observation leads us to the conclusion 
that the scaling form of the PIC persistence probability with nonzero values of the level H is: 

Ppicix, L, 6x, H) = f2{x/L., 5x/L, H/L'^), (15) 

where f2{xi, X2, x^) exhibits a power-law behavior with exponent 0f/c as a function of xi for small xi if a;2 ^ 1 
and — > 0. As the value of X3 is increased, the range of xi values over which the power-law behavior is obtained 
decreases and a more rapid decay of the probability is noticed. 

The predictions concerning the scaling behavior of the spatial persistence probabilities are confirmed by the results 
for the atomistic Kim-Kosterlitz model. The same discussion for Fig. |2 applies to Fig. O where we have shown the 
results for the Kim-Kosterlitz model. We find that 9ss — 0.52 (see Fig. |3Ji), in good agreement with the expected 
value of 1/2, and also that Opic — 0.47, using a rather small simulation with L = 300 and 6x = 1 and sampling 
over the subensemble of sites with height at the average level (see Fig. I^:). As shown in Fig.|3), the SS persistence 
probability obeys the scaling form of Eq. (|14|l . In Fig.^i, we display the results for the measured Ppic for systems 
with different sizes and sampling distances such that Sx/L remains constant and considering two different subsets of 
sampling sites, each subset being characterized by a fixed value of LI/L". These results are in perfect agreement with 
the scaling form of Eq. 11511 . 

Equations (|14|l and (I15() provide a complete scaling description of the SS and FIC persistence probabilities for 
(l+l)-dimensional fiuctuating interfaces belonging to two different universality classes (i.e. EW and KPZ), modeled 
using discrete solid-on-solid models. The associated spatial persistence exponents 9ss and 0f/c are in good agreement 
with the theoretical values 8]. However, these studies do not illustrate the interesting possibility of a dependence of 
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FIG. 3: The spatial persistence probabilities, Pssix) and Pfic[x), for the (l+l)-dimensional Kim-Kosterlitz model which is 
in the KPZ universality class. As in Fig. |21 in panels (a) and (b) we show the data for Pss(x). Panels (c) and (d) display the 
data for Pfic{x). Panel (a): Pssix) for L — 1000, Sx = 1. Panel (b): Finite-size scaling ot Pss{x, L,Sx). Three probability 
curves are obtained for three different sample sizes with the same value for the ratio 5x/L = 1/500. Panel (c): Scaling of 
Pfic(x, L, Sx, H), obtained by sampling over the lattice sites with H = 0, for three different values (same as those in panel (b)) 
of L and Sx. Panel (d): Scaling of Pfic{x, L, Sx, H) for three different sample sizes with the same value for the ratio Sx/L, 
sampling over two subsets of lattice sites with the same value of H/L" {a — 0.5): l/\/300 (upper plot) and 3/\/300 (lower 
plot). 



the persistence exponent on the samphng procedure used in the selection of the initial sites used in the calculation of 
the persistence probability: the two persistence exponents 6ss and Opjc have the same value for (1+1) -dimensional 
EW and KPZ interfaces. We present and discuss below the results for a model where these two exponents have 
different values. 



B. EW equation with colored noise 



In order to measure the spatial persistence and survival probabilities in this system, we have applied the steps 
described above on systems of sizes ~ 2^ — 2^", using 100-400 independent realizations for averages. While the 
calculation of Pgg{x) and Sss{x) involves the same method as the one used in the case of the solid-on-solid models, 
for measuring Pfic{x) and Sfic{x) we have selected the subensemble of lattice sites whose heights h(xj,to) at time 
to ^ satisfy the condition (h) — w/2 < h{xj,to) < (h) + w/2, where (h) is the spatial average of the height at time 
to. The width w of the sampling window has to be chosen to be much smaller than the amplitude of the interface 
fluctuations, but large enough to include a relatively large fraction of the total number of sites in order to ensure 
adequate statistics. Under these circumstances we have computed the fraction of these selected sites which do not 
reach the "original" height h{xj,to) (in the case of persistence probability) or the average height level (h) (in the case 
of survival probability) up to a distance x from the point Xj. The numerical results for these probabilities, along with 
a finite-size scaling analysis of their behavior, are shown in Figs.^and[Sl 

We find that both SS and FIC spatial persistence probabilities for (l-l-l)-dimensional interfaces described by the 
EW equation with colored noise exhibit the expected power-law behavior as a function of x, as shown in Fig. 01 
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FIG. 4: Spatial persistence and survival probabilities for the EW equation with spatially correlated noise. Panel a): Pss{x) 
and Sss{x) using a fixed system size L — 2^, two values of the noise correlation parameter (p — 0.1 and 0.2) and sampling 
distance 6x — 1. Panel b): Ppic{x) and Sfic{x) (inset), using the same parameters as in panel a), and sampling initial sites 
from a band of width w — 0.10 centered at the average height. The straight lines drawn through the data points in these 
double-log plots represent power-law fits. 



while the SS survival probability shows a more complex a;-dependence (see Fig. 0Ji). Further work is needed in order 
to understand the behavior of Sssi^)- When a relatively small system with size L = 2^ is used, the numerical 
results for the spatial persistence exponents extracted from the power-law fits shown in Fig. 0] (for p = 0.1, we obtain 
9ss — 0.43 and Opic — 0.38, while for p = 0.2, the exponent values are found to be 0ss — 0.37 and 9pic — 0.31) 
appear to be affected by finite-size effects. Specifically, the values of 9ss extracted from fits to the numerical data 
are systematically larger than the theoretically expected values, 9ss = 0.4 for p = 0.1 and 0.3 for p = 0.2 (see 
Eq. IH))). Similar deviations from the analytical results are also found for the usual dynamical scaling exponents a 
and (3. We have checked that simulations of larger samples bring the measured values of the exponents closer to the 
expected values, but the convergence is rather slow. These finite-size effects become more pronounced as the noise 
correlation parameter p is increased. In Fig. 0] we show the results for p = 0.1 and p = 0.2, but we have verified from 
simulations with larger values of p that the difference between the expected and measured values of 9ss increases as 
p is increased. This is expected because the spatial correlation of the noise falls off more slowly with distance as p 
is increased, thereby making finite-size effects more pronounced. Another possible source of the discrepancy between 
the numerical and exact results for the exponent 9ss is the spatial discretization used in the numerical work. The 
effects of using a finite discretization scale Ax on the observed scaling behavior of continuum growth equations in the 
steady state have been studied in Ref. p9| where it was found that the effective value of the roughness exponent a 
obtained from calculations of the local width using a finite Ax is smaller than its actual value. Since 9ss = 1 — a, 
the values of 9ss obtained from our calculations with Ax = 1 are expected to be larger than their exact values. Our 
results are consistent with this expectation. As shown in the inset of Fig.^jD, the FIC survival probability Sfic{x) 
behaves similarly to Pfic{x) for both p = 0.1 and 0.2, exhibiting a power-law decay with an exponent (of 0.38 and 
0.33 for p ~ 0.1 and 0.2, respectively) that is very close to 9fic- This is consistent with the expectation that the FIC 
persistence and survival probabilities should become identical as the width parameter w used in the selection of initial 
sites approaches zero (in this limit, both persistence and survival probabilities measure the probability of not crossing 
the average height). Finally, we point out that both SS and FIC exponents obtained from the numerical study exhibit 
the correct trend, increasing in magnitude as p decreases. Also, the measured FIC spatial persistence exponents 
satisfy the constraint 1/4 < 9fic < 1/2. Our numerical results also confirm the interesting theoretical prediction 
that the SS and FIC spatial persistence exponents are different for the EW equation with spatially correlated noise. 

We have found that the scaling forms of Eqs. (|14|l and H15|) also provide a correct description of the numerically 
obtained persistence and survival probabilities for the EW equation with spatially correlated noise. This is illustrated 
in Fig.|31 In Fig. we show that the results for Pss{x, L,5x) obtained for different values of L and 5x fall on the 
same scaling curve when plotted against x/L if the ratio 6x/L is held fixed. This is precisely the behavior predicted 
by Eq. H14|) . As shown in Fig. [5Jd, the data for Pfic{x, L,Sx,w) also exhibit good finite-size scaling collapse if 5x is 
varied in proportion to L and the width w of the sampling band is increased in proportion to L". This is in perfect 
analogy with the scaling behavior of the FIC persistence probability for the discrete stochastic models discussed in 
Sec. IIV Al with the variable w playing the role of H in Eq. (|15|l . This suggests that the scaling behavior of the FIC 
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FIG. 5: Finite-size scaling of the persistence probabilities, Pssix) and Pfic{x), and the FIC survival probability Sfic{x) 
for the EW equation with spatially correlated noise. The noise correlation parameter is p = 0.2 and the sampling interval 
8x takes three different values. Panel a): The SS persistence probability Pss{x, L, 5x) for three different sample sizes with a 
constant ratio 5x/L — 1/2^. Panel b): The FIC persistence probability Pfic{x, L, Sx,w) with fixed values of the quantities 



5x/L (= 1/2**) and w/L"' 
Sfic{x, L, 5x, w). 



0.1/2'^"), where a = 0.7 Inset: Same as in the main figure, but for the FIC survival probability 



persistence probability in the continuum EW equation is of the form 

Ppicix, L,6x,w) = fsix/L, 6x/L, w/L"), (16) 

where the function /s has the same characteristics as /2 in Eq. (|15|l . A similar scaling description also applies to 
Sfic{x), as shown in the inset of Fig. 03. This scaling description should be useful in the analysis of experimental 
data on equilibrium step fluctuations JiJ because the images obtained in experiments provide the values of a real 
"height" variable (position of a step-edge) at discrete intervals of a finite sampling distance 5x. 



V. SUMMARY AND CONCLUDING REMARKS 



In this study, we have analyzed the spatial first-passage statistics of fiuctuating interfaces using the concepts 
of spatial persistence and survival probabilities. Specifically, we have presented the results of detailed numerical 
measurements of the SS and FIC spatial persistence probabilities for several models of interface fluctuations. Results 
for the spatial survival probabilities are also reported. These results confirm that the concepts of persistence and 
survival are useful in analyzing the spatial structure of fluctuating interfaces. The exponents associated with the 
power-law decay of the spatial persistence probabilities as a function of distance x are valuable indicators of the 
universality class of the stochastic processes that describe the dynamics of surface fluctuations. Our results for these 
exponents for (l-l-l)-dimensional interfaces in the EW and KPZ universality classes are in good agreement with the 
corresponding analytic predictions. We have also obtained analytic results for the spatial persistence exponents in the 
(l-l-l)-dimensional EW equation with spatially correlated noise, and reported the results of a numerical calculation 
of the persistence and survival probabilities in this system. While the numerical results show strong finite-size effects, 
the qualitative trends predicted by the analytic treatment are confirmed in the numerical work. In particular, the 
numerical results show evidence for an interesting theoretically predicted difference between the persistence exponents 
obtained for two different ways of sampling the initial points used in the measurement of the spatial persistence 
probability. We also find that the steady-state survival probability has a complex spatial behavior that requires 
further investigations. In the past, there has been some confusion in the literature about the distinction between the 
persistence and survival probabilities |l5j |. Our study shows that these two quantities are very different in the SS 
situation, whereas the distinction between them essentially disappears in the FIC situation. 

The numerical results reported here are for models that exhibit "normal" scaling behavior with the same local 
and global scaling properties of interface fluctuations. There are other models of interface growth and fluctuations 
that exhibit "anomalous" scaling j23|, for which the global and local scaling properties are different. In such models, 
the "global" roughness exponent ag that describes the dependence of the interface width in the steady state on the 
sample size L {W{tQ,L) oc L"" for to S> L^) is different from the "local" exponent ai that describes the x-dependence 
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of the height-difference correlation function g{x) = {[h{x + Xq^Iq) — /i(a;o: ^o)]^)^^^ in the steady state [t^ ^ L^) for 
small X {g{x) oc a;"' for x <^ L). The exponent ag is greater than unity (the steady-state interface is "super- rough" ) 
in such cases, whereas the local exponent ai is always less than or equal to unity. It is interesting to enquire about 
the behavior of the spatial persistence probabilities in such models. The numerical results reported in the preceding 
sections show that the steady-state persistence probability Pssix) exhibits a power-law decay in x only for values 
of X that are much smaller than the sample size L. Since the roughness of the steady-state interface of super-rough 
models at length scales much smaller than L is described by the local exponent ai , we expect the steady-state spatial 
persistence probability in such models to exhibit a power-law decay with exponent 9ss = 1 — 0; for x <^ L. For 
example, the one-dimensional MuUins-Herring model is super-rough with ag = 3/2 and a; = 1. For this model, 
the above argument suggests that the steady-state spatial persistence exponent 9ss is equal to 0, which agrees with 
the exact result reported in Ref. 

An important feature of our investigation is the development of a scaling description of the effects of a finite system 
size and a finite sampling distance on the measured persistence probabilities. We have also shown that the dependence 
of the FIC persistence and survival probabilities on the reference level H (in atomistic models) or the width w of the 
band (in continuum models) used in the selection of the subset of sampling sites is described by a scaling form. These 
scaling descriptions would be useful in the analysis of experimental and numerical data on fluctuations in spatially 
extended stochastic systems. 

Some of the numerical results reported here (such as the behavior of the SS survival probability and the forms of the 
scaling functions that describe the dependence of the persistence probabilities on the parameters L, Sx and H or w) 
should be amenable to analytic treatment, especially for the EW equation with white noise, whose spatial properties 
can be mapped to the temporal properties of the well-known random walk problem. Further work along these lines 
would be very interesting. The spatial persistence and survival probabilities considered here should be measurable in 
imaging experiments on step fluctuations |^|^. Such experimental investigations would be most welcome. 
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